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Abstract 

Driven by the loss of energy, isolated rotating neutron stars (pulsars) are gradually slowing down to lower 
frequencies, which increases the tremendous compression of the matter inside of them. This increase in 
compression changes both the global properties of rotating neutron stars as well as their hadronic core 
compositions. Both effects may register themselves observationally in the thermal evolution of such stars, 
as demonstrated in this Letter. The rotation-driven particle process which we consider here is the direct 
Urea (DU) process, which is known to become operative in neutron stars if the number of protons in the 
stellar core exceeds a critical limit of around 11% to 15%. We find that neutron stars spinning down from 
moderately high rotation rates of a few hundred Hertz may be creating just the right conditions where 
the DU process becomes operative, leading to an observable effect (enhanced cooling) in the temperature 
evolution of such neutron stars. As it turns out, the rotation-driven DU process could explain the unusual 
temperature evolution observed for the neutron star in Cas A, provided the mass of this neutron star lies 
in the range of 1.5 to 1.9 M Q and its rotational frequency at birth was between 40 (400 Hz) and 70% (800 
Hz) of the Kepler (mass shedding) frequency, respectively. 
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Numerical cooling simulations of neutron stars 
allow one to probe the inner structure of these 
objects, and the properties of ultra-dense matter 
[113, i, i, i, H, 0, II M ■ In the standard treatment, 
neutron star cooling calculations are carried out for 
stellar core compositions which are frozen-in, that 
is, compositions which do not change with time, as 
it is the case for non-rotating neutron stars. The 
situation may be very different for isolated rotating 
neutron stars, which are spinning down due to mag- 
netic braking. Such stars can experience drastic 
density changes during spin-down, causing the for- 
mation of novel states of matter in their cores. Ex- 
amples of which are the formation of a mixed phase 
of quarks and hadrons, of pure quark matter, or a 
condensate made of bosons 2, Ell El [IS El- An- 
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other intriguing possibility concerns the rotation- 
driven changes in the number densities of neutrons, 
protons and leptons in the cores of rotating neutron 
stars, which can have important observable conse- 
quences for the thermal evolution of such objects, 
as shown in this Letter. 

Computing the thermal evolution of rotating 
neutron stars is considerably more complicated 
than computing the cooling of non-rotating neutron 
stars. The reasons are several. First, stellar rota- 
tion requires solving Einstein's field equations for 
rotationally deformed fluid distributions [13, 11|. 



which renders the problem 2-dimensional. Second, 
the general relativistic frame dragging effect im- 
poses an additional self-consistency condition on 
Einstein's field equations. Third, an extra self- 
consistency condition is encountered when calculat- 
ing the general relativistic Kepler (mass shedding) 
frequency of a rotating neutron star. Fourth, the 
thermal transport equations need to be solved for 
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Figure 1: (Color online) Critical temperature of supcrfluid 
neutrons in 1 Sq and 3 P2 states as a function of baryon num- 
ber density, (no denotes the density of ordinary nuclear mat- 
ter.) 



general relativistic, non-spherical fluids that may 
experience anisotropic heat transport. Our study 
accounts for the first three effects. The fourth ef- 
fect, the anisotropic transport of heat inside neu- 
tron stars, is reserved for a separate future study 

Ell- 

In this Letter, we consider the cooling behav- 
ior of isolated rotating neutron stars. The num- 
ber of baryons (i.e., the so-called baryon mass) of 
such stars remains unchanged during spin-down. 
As known from earlier studies 1^, IH 12 [ > the spin- 
down driven gravitational compression can cause 
substantial changes in the core compositions of such 
neutron stars. A striking example of this is the di- 
rect Urea (DU) process, which describes the direct 
(no bystander particle necessary) transformation 
of neutrons to protons, electrons and anti-electron 
neutrinos in the cores of neutron stars, according 
to 



n-> p + e +u e , (1) 

Reaction (JTJ) can only occur if the proton fraction 
exceeds between 11% to 15% [lq |. The neutrino 
luminosity associated with this reaction dominates 
over those of other neutrino emitting processes in 
the core [H, [l8[ . One may thus speculate that once 
([1} becomes operative in rotating neutron stars, it 
could speed up the cooling of such stars substan- 
tially, leading to an observational signature in the 
thermal evolution of such objects. The pairing of 
neutrons in the 1 So singlet and 3 P2 triplet states, 
which reduces the neutrino emission from neutron 
stars, would not change this behavior, as will be dis- 



cussed later in this Letter. Following the approach 
of 0, [20[, the critical temperatures assumed for 
these states are shown in Fig. [T] 

Considering the conservation of momentum and 
the superfluidity suppression, we show in Fig. [5] the 
intensity of the DU process (defined here as the re- 
duction factor introduced by pairing, i.e. "1" means 
no reduction, "0" total suppression) in the core of 
rotating neutron stars for different stellar masses, 
M, and rotational frequencies, f2. One sees that, 
depending on frequency (and thus central density) 
and temperature (which affects the level at which 
the DU process is suppressed), the DU process oper- 
ates at very different intensities. Hence, as a rotat- 
ing neutron star evolves from a given "initial" state 
(given frequency, mass, and temperature) in time, 
the intensity of the DU process at its core can vary 
substantially. Also shown in Fig. [5] are several evo- 
lutionary tracks of neutron stars which spin down 
because of the emission of magnetic dipole radia- 
tion. The intensity of the DU process of stars along 
these tracks can vary significantly too, depending 
on the values of the star's initial parameters, which 
ought to affect their thermal evolution. 

In passing we note that the rotationally driven 
compositional changes in the cores of rotating neu- 
tron stars may also occur, in reverse however, in the 
cores of neutron stars in binaries, which are spun 
up by mass accretion. The treatment of such neu- 
tron stars is more complicated, however, since both 
the frequency as well as the star's mass are chang- 
ing. In that case the evolutionary tracks shown in 
Fig. [5] will follow a different path, but the overall 
modification of the intensity of the DU process may 
be expected to hold. 

The results shown in Fig. [2] are generic for any 
EoS which predicts a sufficiently large number of 
protons in the cores of neutron stars so that re- 
action ((TJ) becomes possible. The EoS used in 
this Letter is a parametrized version of the Akmal- 
Pandharipande-Ravenhall (APR) EoS [HIH] ( and 
references therein). For this parametrization, the 
binding energy consists of a compression term and 
a symmetry energy term, 



BE = E u 



2-5 



1 



+ S u 1 (1 - 2x p ) 2 , (2) 



where Eq is the saturation energy, So the symmetry 
energy at saturation, u = n/no with n the baryon 
number density and no the baryon number density 
at saturation, and x p the proton fraction. The pa- 
rameters 5 and 7 can be used to control the nuclear 
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Figure 2: (Color online) Intensity of the direct Urea (DU) process in rotating (frequency f2) neutron stars, for different stellar 
temperatures, T. M denotes the gravitational mass. The intensity ranges from (plain white areas) to 100% (red). No stars 
are allowed above the curve labeled Kepler limit because of stellar mass shedding. The vertical curves show evolutionary tracks 
of isolated rotating neutron stars, whose baryon mass (Mo) during spin-down remains constant. One sees that for several 
stars evolving along these tracks the DU process (i.e., reaction |[TJ) becomes possible at lower rotation rates. The color coding 
indicates the intensity of the DU process. 



incompressibility and the symmetry energy, respec- 
tively. The maximum mass of a neutron star is thus 
sensitive to changes in S, and the threshold den- 
sity for the onset of the DU process is sensitive to 
changes in 7. Details of the parametrization of the 
APR EoS have been discussed in [22j. The cooling 
of neutron stars described by this EoS were first ex- 
plored in ^23]. Very recently, this EoS was used to 
put constraints on the properties (superfluid gaps) 
of ultra-dense matter using the temperature data 
observed for the neutron star in Cas A [1, H, 24 1 . 

Next, we explore the impact of a varying pro- 
ton fraction, driven by stellar rotation, on the ther- 
mal evolution of isolated rotating neutron stars. In 
this study we consider the following cooling pro- 
cesses: the direct Urea, the modified Urea, and 
Bremsstrahlung processes in the stellar core. For 
the crust, we consider Bremsstrahlung, e + e~ pair 
annihilation, and plasmon decay processes. A com- 



prehensive overview of the neutrino emission pro 



cesses in neutron stars can be found in [101 
The boundary conditions are determined by the lu- 
minosity at the stellar center, 1(0) = 0, and at the 
surface l(R) = L s [26|,[27|. To establish a connec- 
tion between the rotational frequency of neutron 
stars and their respective cooling stages, we model 
the spin-down rate according to 



n = - k n r 



(3) 



where f2 denotes the neutron star spin frequency, 
and K and n are constants The constant K 

can be parametrized as (for n = 3) K = 1.55 x 
lCT 17 x Rf Bl 2 /(M/M )s, where B 12 is the sur- 
face magnetic field in units of 10 12 G, Rio the ra- 
dius in units of 10 km, and a the inclination an- 
gle of B (note that the moment of inertia was ap- 
proximated by 2MR 2 /h). Using the frequencies 
and time derivatives of observed pulsars available 
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at the ATNF pulsar data base [28|, we find that 
10 -21 < K < 10~ n (in the same units as above, 
i.e., s) for a braking index n = 3. The latter value 
corresponds to neutron stars which spin-down due 
to the emission of magnetic dipole radiation, as con- 
sidered in this Letter. By integrating Eq. ([3]) one 
obtains the frequency as a function of time. With 
the aid of Eq. ([3]), we can now trace the thermal 
evolution of isolated, spinning-down neutron stars. 

The results are shown in Figs. [3] and |U It is 
evident from Fig. [3] that the temperature drops 
dramatically when the DU process ([T]) becomes 
operative during spin-down. The characteristic 
time scale of the thermal coupling between the 
core and the crust of non-rotating stars is tc — 
(Ai?/lkm) 2 (l - 5.04 x 10~ 2 (M/M Q )/R w )- 3 / 2 h, 
where AR is the crust thickness, and t\ is a con- 
stant with dimension of time, that depends on the 
microscopic properties of the star, and whose value 
is ~ 30 years [29j. For a traditional neutron star 
tc = 50 — 150 years [29j], and in this work we set 
tc = 100 years . In the case of spinning-down neu- 
tron stars considered here, we must also consider 
the relaxation time of the spin-evolution, which we 
define as the time it takes for the star to spin down 
to half of its original frequency. This time scale 
is given by ts — 3/(2KHq), where fio denotes the 
birth frequency of the object. The parameter 







ts_ _ M/M & (1 - 5.04 x 10~ 2 (M/M Q )/R 10 ) 3/2 

TC 



1.03 x 10- 17 Rj B 2 2 sin 2 aAR 2 



(4) 

can then be used to determine the impact of the 
spin-down compression on the cooling of the star. 
Evidently this analysis is only valid for objects that 
are born in the white regions of the figures shown 
in @, and thus may cross-over into the DU region 
(red shaded areas) during their evolution. For the 
EoS studied here, this is limited to objects whose 
non-rotating gravitational masses are > 1.5 M e . 
The following cases emerge: I. If /3 < 1, the star 
spins down to low frequencies (fl — > 0) before the 
core and crust are thermally coupled. This causes 
the cooling of the star to be similar to the cooling 
of a spherically symmetric star of same mass. II. If 
P ~ 1, the star is still rotating at relatively high fre- 
quencies (~ rio/2) when the core and crust become 
thermally coupled. The thicker crust of such a star 
allows the core and crust to couple more quickly, 
effectively speeding up the cooling. III. If /? > 1, 
the star keeps rotating at high frequencies for a long 
time, which delays the onset of the DU process (JTJ) - 



We note that, as was the case for (3 ~ 1, the star 
also becomes isothermal at an earlier time. How- 
ever the late onset of the DU process decouples the 
core and crust one more time, which yields to a 
second thermal coupling, characterized by a sharp 
drop in temperature at a later time. IV. If (3 — > oo, 
the spin-down relaxation time is much greater than 
the time scale of the core-crust thermal coupling. 
This means that during most of the cooling period, 
the object will remain in the high frequency do- 
main, and the onset of the DU process may never 
be achieved, leading to a slow stellar cooling. The 
different scenarios I through IV are graphically il- 
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Figure 3: (Color online) Neutron star temperature at infinity 
as a function of stellar age, for a sample neutron star with 
baryon mass of 1.5 Mq. (The gravitational mass at zero 
rotation is 1.4 Mq.) Depending on the value of K, four dis- 
tinctly different cooling epochs emerge, which are discussed 
in the text. 

lustrated in Fig. [31 The curves underline the impor- 
tant role of /3 for the classification of the cooling be- 
havior of spinning-down neutron star. If superfluid- 
ity effects were considered, one should expect that 
the temperature reductions for /3 ~ 1 and /3 > 1 
are less pronounced. The actual magnitude of the 
reduction will depend on the microscopic pairing 
model used. 

We conclude this Letter with a discussion of the 
temperature data observed for the neutron star in 
Cas A over a 10 year period [30j. The rapid cool- 
ing of this neutron star has been explained recently 
through the presence of superfluidity in dense mat- 
ter [8j, |{|. Our study indicates that the tempera- 
ture evolution of this neutron star could also be ex- 
plained by the late onset of the DU process, if one 
assumes that the neutron star was created with an 
initial rotational frequency somewhere between 40 
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and 70% of its mass shedding frequency (see Fig. 
2]), depending on the unknown mass of this neutron 
star. Masses between 1.2 and 2.1 M Q (99% confi- 



2.12x10 
2.10xlO f 
2.08x10* 
2.06X10 6 
2.04xlO f 
2.02X10 6 
2.00xlO f 



A K = 6.90sl0 .AM = 5.5x10 ,n /£l K = 0.67 

\ K = 3.34 x 10 ' 6 . AM = 6.0 x 10"' , n o /n K = 0.42 

♦ — ■ K=1.84xl0"\AM=1.0xl0" B , 0^ = 0.37 

\ \ M = 1.5 M , B = 21.08 (SF) 

\ \ 




310 



320 330 340 

Age (years) 



350 



Figure 4: (Color online) Cooling simulations reproducing the 
temperatures observed for the neutron star in Cas A over 
a ten-year period. K is defined in the text, AM denotes 
the mass (in unit of the mass of the neutron star) of the 
accreted envelope, and f!o/^it denotes the star's birth fre- 
quency relative to the Kepler frequency. Superfluidity and 
Pair-Breaking Formation of nucleons is taken into account 
in the curve labeled SF. 

dence) have been established for this neutron star 
The simulations shown in Fig. H] have thus been 
performed for two sample masses, 1.5 and 1.9 Mq. 
The mass of the accreted envelope, due to the fall- 
back of matter after the supernova explosion 27j, is 
assumed to be AM/M = 5.5 x 10 _n and 6.0 x 10" 9 , 
in accordance with (3l| . Our calculations predict a 
birth frequency (Oq) for the 1.5 Mq star of 42% 
(f2o = 400 Hz) of the Kepler frequency if the stars 
is made of non-superfluid matter. This figure drops 
down to 37% if neutron superfluidity in the 1 Sn 
and 3 P2 channels and Pair-Breaking Formation [8j 
is included. It is an open issue whether proton su- 
perconductivity takes place in the core of a neutron 
star. If it does, it may affect the thermal evolution 
of the object [8j , depending on the magnitude and 
density dependence of the superfluid gap and its 
critical temperature. In the case that proton su- 
perconductivity extends to very high densities, and 
has a high critical temperature, it might suppress 
the impact of the DU process on the thermal evo- 
lution of a neutron star. 

We note that on the basis of theoretical calcula- 
tions of the mapping between initial and final neu- 
tron star spins, likely spin-down mechanisms, and 
observational constraints from pulsars and super- 
nova energetics, it has been argued |32j that the ini- 



tial pre-collapse central iron core periods inside of 
massive stars, which give birth to neutron stars and 
pulsars, are on average greater than ~ 50 seconds. 
The associated neutron stars would then be born 
with rotation periods greater than around 10 ms 
[32| . This fi gure agrees with the birth periods es- 
tablished in 33j. If the central iron core periods 



should be smaller by a factor of 5, the birth periods 
of neutron stars would drop down to a millisecond. 
In our work, initial birth periods of 2.5 ms (400 Hz) 
and 1.25 ms (800 Hz) are required to fit the thermal 
data of the neutron star in Cas A. Such rotational 
periods are on the small end of probabilities, but 
given the many poorly understood issues (estimated 
iron core spin rates, angular momentum profiles, 
progenitor mass, general relativistic effects), which 
complicate the study of stellar evolution with rota- 
tion, the birth periods considered in our study are 
in the realms of possibility. 

A further aspect concerns the energy loss dur- 
ing spin-down for the assumed birth frequencies. 
For our model scenario, we find that E = I tin ~ 
10 39 erg/s. This figure would imply that the neu- 
tron star in Cas A ought to emit a pulsar wind, 
which, however, has not been observed. If the pul- 
sar wind should not exist, the birth frequencies of 
the neutron star in Cas A may have been some- 
what smaller than the values assumed in this paper 
for neutron stars made of non-superfluid matter. 
Interestingly, the inclusion of superfluidity among 
neutrons and Pair-Breaking Formation reduces the 
region in the star where the DU process takes place, 
shifting the onset of the DU process toward lower 
rotational frequencies (see curve labeled "SF" in 
Fig. H]). 

The results shown in Fig. [4] were obtained by as- 
suming value of K in the range of ~ 10~ 16 and 
~ 10~ 18 . The value of K is connected to the mag- 
netic field through K oc Rf B% 2 /(M/M e ). It is 
known from observations that the upper limit for 
the magnetic field of the neutron star in Cas A is 
2 x 10 11 G, which would imply K ~ 10~ 18 (assum- 
ing a stellar mass of M =2.0 Mq and a radius of 
R = 15 km), in agreement with one of the parame- 
ter combinations considered in Fig. 0] (blue curve) . 

We also note that given the different nature of 
the processes employed to explain the temperature 
drop (late onset of the DU process in our case, 
and pair breaking/formation in references [H, Q), 
continuous observation of the neutron star in Cas 
A may decide which model is favorable. A con- 
tinuous analysis of the slope of the temperature 
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(s = d log 10 Ts/d log 10 1) might allow one to do 
that. While the slope found by us (s — —1.5) should 
remain at this value while the surface temperature 
drops drastically, in the model described in [8], it 
should go to the asymptotic value found by the au- 
thors (s = —1/12), before the surface temperature 
changes drastically. 

Our results show that rotation-driven repopu- 
lation processes can be of very great importance 
for the thermal evolution of spinning-down neu- 
tron stars. The reason for this is two-fold: first, 
the geometry of the object is modified as they spin 
down, which changes the stellar surface gravity and 
therefore the stellar surface temperature [27]. Sec- 
ond, the spin-down renders neutron stars gradually 
more dense, which changes the stellar core composi- 
tions. We also discovered that the ratio of the spin- 
down relaxation time to the core-crust coupling 
time, /3, emerges as a most valuable parameter that 
serves to classify the cooling behavior of rotating 
pulsars. Furthermore, the methodology discussed 
in this Letter can be generalized to more elabo- 
rated scenarios like accreting (X-ray) neutron stars, 
and possibly objects undergoing more rapid braking 
mechanisms like r-modes. Evidently, in such cases, 
re- heating processes are important, although they 
do not invalidate the conclusions drawn by us in 
this Letter, i.e. that particle repopulation plays an 
important role for the thermal evolution of neutron 
stars. 

In summary, the key point put forth in this let- 
ter is that the particle composition in the core of 
an isolated rotating neutron star changes due to the 
spin-down compression of the star caused by mag- 
netic braking. This compression is a robust physical 
phenomenon, experienced by any isolated rotating 
neutron star. This scenario is applied to the neu- 
tron star in Cas A. We show that the remarkable 
temperature drop of this neutron star can be ex- 
plained by the onset of the rotation-driven direct 
Urea process. We also show that continuous obser- 
vation of the thermal evolution of this neutron star 
over the next ten to twenty years will allow us to 
determine whether our explanation of the fast cool- 
ing of this neutron star is correct, or whether other 
physical processes (as suggested in 0, Q) are re- 
sponsible for the star's dramatic temperature drop. 
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